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When visualised as an operation on the Bloch sphere, the qubit "pi-over-eight" gate corresponds to 
one-eighth of a complete rotation about the vertical axis. This simple gate often plays an important 
role in quantum information theory, typically in situations for which Pauli and Clifford gates are 
insufficient. Most notably, when it supplements the set of Clifford gates then universal quantum 
computation can be achieved. The "pi-over-eight" gate is the simplest example of an operation from 
the third level of the Clifford hierarchy (i.e., it maps Pauli operations to Clifford operations under 
conjugation). Here we derive explicit expressions for all qudit (d-level, where d is prime) versions 
of this gate and analyze the resulting group structure that is generated by these diagonal gates. 
This group structure differs depending on whether the dimensionality of the qudit is two, three 
or greater than three. We then discuss the geometrical relationship of these gates (and associated 
states) with respect to Clifford gates and stabilizer states. We present evidence that these gates are 
maximally robust to depolarizing and phase damping noise, in complete analogy with the qubit case. 
Motivated by this and other similarities we conjecture that these gates could be useful for the task 
of qudit magic-state distillation and, by extension, fault-tolerant quantum computing. Very recent, 
independent work by Campbell, Anwar and Browne confirms the correctness of this intuition, and 
we build upon their work to characterize noise regimes for which noisy implementations of these 
gates can (or provably cannot) supplement Clifford gates to enable universal quantum computation. 



I. BACKGROUND AND MOTIVATION 

The qubit "pi-over-eight" gate, 17^/8, 

v.» - ( '7 A ) w 

plays a special role in a number of quantum informational 
tasks. The Gottesman-Knill theorem jlj tells us that a 
circuit using only Clifford gates and Pauli measurements 
(i.e., a stabilizer circuit) is insufficient for universal quan- 
tum computation (UQC). Whilst technically, adding the 
ability to perform any single-qubit non-Clifford gate is 
sufficient for obtaining UQC, one typically sees that the 
is chosen as the most natural and easiest with which 
to work [2]. In measurement-based scenarios, supple- 
menting Pauli measurement directions with an additional 
rotated (by U v /s) measurement basis can also enable the 
performance of new tasks. For example, performing Pauli 
measurements on the Bell state (|00) + |ll))/\/2, or any 
two-qubit stabilizer state, does not exhibit any better- 
than-classical performance in a nonlocal CHSH game 
[31 H] , whereas the introduction of the new rotated ba- 
sis enables the optimum quantum advantage. Measure- 
ments in the aforementioned rotated basis appear to arise 
naturally in other quantum- informational tasks too e.g., 
universal blind quantum computation [5], where Pauli 
measurements and operators would be insufficient. 
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Arguably, much of the utility of this gate arises from 
its close relationship with the Clifford group, whilst still 
not being a member of the group. In fact, U^/g is the sim- 
plest meaningful example of a gate from the third level 
of the Clifford hierarchy (defined later), the first two lev- 
els of which correspond to Pauli gates and Clifford gates. 
Operations from the Clifford hierarchy have properties 
that make them suitable for teleportation-based UQC 
[B], transversal implementation (see below), learning an 
unknown gate [7j, or secure assisted quantum computa- 
tion [8] . Gates from higher levels of the Clifford hierarchy 
are also related to so-called semi-Clifford gates [91411). 

It is now known |12j that no quantum error-correcting 
code allows for transversal (i.e., bitwise and manifestly 
fault-tolerant) implementation of a universal set of gates. 
Typically, many stabilizer codes (e.g., CSS codes) enable 
transversal implementation of the complete set of Clifford 
gates, but only these gates. An intriguing exception is 
provided by the [[15,1,3]] punctured Reed-Muller code 
that allows for transversal implementation of the U n /g 
gate, but not the complete set of Clifford gates. Zeng 
et al. |13) provide quite a detailed discussion on the 
relationship between Reed-Muller codes and transversal 
gates from the third (or higher) level of the Clifford hi- 
erarchy. It appeared that the transversal U w /s property 
of the [[15,1,3]] code was useful in the derivation of a 
magic state distillation (MSD) routine by Bravyi and Ki- 
taev |14) , a routine that iteratively distills increasing pure 
copies of the state \ipu w / a ) — ^7r/sl+) ^ |0) + e'*|l). 

Independently of the foregoing discussion, it is also no- 
table that U-x/8 is remarkable due its geometrical relation- 
ship with the set of Clifford gates - it is the single-qubit 
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unitary U € SU(2) that is farthest outside the convex 
hull of Clifford operations. The associated single-qubit 
state \ipu„/ a )i mentioned previously, is also remarkable 
in its convex-geometrical relationship with Pauli eigen- 
states. Furthermore, these geometrical relationships have 
ramifications for the amount of noise that can be toler- 
ated by imperfect implementation of U^/s or imperfect 
preparation of \tpu^/ s }-> a scenario that arises naturally in 
any fault-tolerant UQC proposal that uses magic state 
distillation. 

The state \ipu„ /s ) oc |0) + e 4 T|l) is already known in 
quantum information theory as \H) - a qubit magic state 
as defined in [21 [15]. In addition, the most non-stabilizer 
qubit state \T) [Mj, in a convex-geometrical sense, is also 
a magic state. Moreover, both \H) and \T) are eigenvec- 
tors of Clifford gates. The importance of geometrically 
significant states and gates to qubit-based fault-tolerant 
UQC provided the motivation in [T5] to find the most 
robust qudit states (for all prime dimension) and qudit 
gates (for p £ {2, 3, 5, 7}). The maximally robust states 
(analogous to \T)) were also found to be eigenvectors of 
Clifford gates, whereas the maximally robust gates had 
a strikingly simple form, which prompted the question 
whether these gates were related to generalized versions 
of the C/jr/s gate. Here, starting from the condition that 
these gates must be diagonal elements of C 3 , we derive 
generalized versions of LL/g, which we call U v and show 
that these are identical, up to an unimportant factor of a 
Clifford gate, to the maximally robust gates found in [16]. 
We also show that the associated states \i/>u v ) — U v \+) 
are eigenvectors of Clifford gates, and that they obey a 
similar relationship with respect to stabilizer states as 
\H) does. 

As we completed this work we became aware of very 
recent results by Campbell, Anwar and Browne [T7] . 
There, the authors prove the existence of magic state 
distillation protocols (MSD) for all (prime) qudit sys- 
tems, wherein the non-stabilizer states that are distilled 
are states that we have called \ipu v ) here. Moreover, they 
show that U v have a transversal implementation in a fam- 
ily of qudit Reed-Muller codes and they explain why this 
property is useful in MSD. It is hoped that the results 
presented here might aid in the analysis of such qudit 
MSD protocols. More generally, it seems likely that the 
gates U v will find application in other areas of quantum 
information theory, particularly in qudit generalizations 
of qubit-based tasks for which U v / 8 is known to be help- 
ful. 

We begin Sec. [TT] by deriving explicit expressions for 
all qudit U v , and then proceed to analyze the resulting 
group structure where we note an interesting difference 
depending on whether p = 2, p = 3 or p > 3. In Sec. |III| 
we discuss geometrical features of these gates U v and 
associated qudit states \ipu v ) with a particular eye to ap- 
plications in quantum computation. We conclude in sec- 
tion Sec. IIVI with some observations on noise thresholds 
for qudit-based quantum computation. 



II. BASIC MATHEMATICAL STRUCTURE 

A. Generalized Pauli and Clifford Groups 

Throughout, we always assume the dimension p, of a 
single particle, to be a prime number. Generalized ver- 
sions of the familiar cr x and o z Pauli operators, are de- 
fined [18 for p > 2 as 

X\j) = \j + lmodp) Z\j)=J\j) (2) 

where u = e 27Tl / p is a primitive p-th root of unity such 
that XZ = lj~ 1 ZX. In general, products of these Pauli 
operators are often called displacement operators, 

D {x \ z) = t xz X x Z z t = e< p+1 ) li / p = uj 2 ' 1 (3) 

where the format of the subscript (x\z) reminds us of 
the symplectic form (often used in calculations involving 
Pauli operators e.g., error-correcting codes). The Weyl- 
Heisenberg group (or generalized Pauli group) for a single 
qudit is given by 

Q = {t c Dz\x e^,ce Z p ) (Z p = {0, l, . . . , P - l}) , 

(4) 

where x is a 2- vector with elements from Z p ,so that \Q\ = 
p 2 in situations where global phases can be ignored. The 
set of unitary operators that map the Pauli group onto 
itself under conjugation is called the Clifford group, C, 

C = {C eu(p)\cg<y =G}. 

The number of distinct Clifford gates for a single qudit 
system (ignoring global phases) is \C\ — p 3 (p 2 — 1). 

Gottcsman and Chuang [6] introduced the so-called 
Clifford hierarchy, a recursively defined set of gates given 
by 

C k+1 = {U\UdU^ CC fe } (5) 

where C\ is the Pauli group. One obtains nested sets 
of operators, the first two sets of which correspond to 
elements of the Pauli and Clifford groups respectively, 
i.e., 

Q C C C C 3 C . . . (6) 

or equivalently in their notation 

Ci C C 2 C C 3 C . . . (7) 

It is known JS[ [HI [TU] that C3 (and above) does not form a 
group, although the diagonal subset of C3, that we study 
here, does. 

The complete set of Clifford unitaries C C U(p) is cov- 
ered by varying over all F e SL(2, Z p ) and \ £ Z p , 

C = {C {m) I F € SL(2,Z p ),x£% 2 p }, (8) 
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where SL(2,Z p ) is the group whose elements are 2x2 
matrices with unit determinant and matrix elements from 
Z p . The explicit recipe [TS] for constructing a Clifford 



unitary with F = 



.7 5 )i X 



is given by 



C, 



(F\x) 



D( x \ z )Vf 



(9) 



_L 



Ep-l 
fe=0 



Note that when F 



VpD 



ak)(k\ 

I2 we have Vf 



= 0. 



Also 



Vl = D 



FJ->(x\z) v F 
C {Fi\xi) C {F 2 \x-2\ « C, 



(aa;+/32|7a:+5z) 



(^1^1^1+^1x2) 



(10) 

(ii) 



where the proportionality symbol denotes equality mod- 
ulo a global phase. 

The particular case p = of Eq. ^ turns out to be 
particularly relevant to our investigation, and so we note 
that 



det 



y T a~fk 2 
\k=0 



T ^(2P-1)(P-1)P 



= 1 Vp>3, 
= t 2 " 7 p = 3, 



(12) 

(13) 
(14) 



which has ramifications for results in the next section. 
(This peculiarity, for p = 3, was also noted by Zhu [2U] 
in the course of an investigation into Weyl-Heisenberg 
covariant SIC-POVMs.) In particular we will use 



C 



det C 



([$?]|[:]) 



e SU(p) Vp > 3, 

= r 27 for p = 3, 



(15) 



which eventually leads to an unusual group structure for 
qutrit generalizations of U_/g. 

B. Explicit form of qudit gates analogous to U n /% 

For simplicity and, in analogy with the qubit 11^/%- 
gate, we chose U v (our putative higher-dimensional gen- 
eralizations of C/tt/s) to be diagonal in the computational 
basis, so that UD^)W = -D(o|i)- We claim that, for 
p > 3, U v can be written in the following form: 



p-i 

U v = U(v ,v l7 ...) 

j=0 



K e Zp), (16) 



where lu = e p , as usual. Note that det(U v ) — oj^k=o Vk 



so that U v € SU(p) if J2k=ayk — (mod p). Straight 



forward application of Eq. |2]), ^ and ( fl6| gives 

U v D (x]z) Ul = D (xlz) Y / ^ {vk+ ^ Vk) \k)(k\ (17) 



If U v is to be a member of C3 we require the right hand 
side of (17) to be a Clifford gate. Since U v D^\ z ^Ul = 
D(o|z), trivially, we focus on the case U v D(i\ )Ul in order 
to derive explicit expressions for U v . 
Define 7', z' , e' £ Z„ such that 



U v D {m W v =w € ^[i 0]|[i]) ( 18 ) 



The fact that the right hand side of the Eq. ( 18 ) is the 
most general form can be seen by reference to Eq. ^ 
and fl7| ), and also by noting that U € SU(p) implies 
w fe £/ G SU(p), for any integer k. 



Note that the right hand side of Eq. ( 18 ) represents a 
Pauli operator if and only if 7' = 0. Consequently, U v 
must, by definition, be a (diagonal) Clifford operation in 
those cases when 7' = 0. 



To solve the matrix equation Eq. (18), begin by sub- 
stituting Eq. (17) so that 



^dio) E t» lVh+i - Vh) \w\ =Lj£ ' c ([\° 1 ]\[\]) { 1 " ] 



and use Eq. ^ to obtain 



D 



p-i 



(20) 



fc=0 



After canceling common factors of -Dmo) one is left with 
an identity between two diagonal matrices, so that 

u vh+i-vi. = ujt' T z ' uj kz ' T k2 ~<' (Vfc E Z„), (21) 



or equivalently, using r = uj 2 



v k+ \ - v k = e' + 2~ l z' + kz' + 2- 1 k 2 j'. 



(22) 
(23) 



This gives the recurrence relation 

= v k + k(2- 1 kj l + z') + 2~V + e'. 
With a boundary condition vq = 0, we can solve to obtain 

Vk = 7^(7' + H6z! + (2k - 3)7')) + ke' (24) 

where factors like 12 _1 are understood to be evaluated 
modulo p. 

For example, with p = 5 and choosing z' = 1,7' = 4 
and e' = 0, we get 



v = (v 0l v 1 ,v 2 ,v 3 ,v i ) = (0,3,4,2, 1) 
>U v (z',i,e') 



/I 














\ 





47ri 

e 5 



















2tvi 

e b 



















47Ti 

e 5 







Vo 











27Ti 

e 5 


/ 



(25) 
(26) 
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It can be shown that the powers of oj along the diagonal 
of U v sum to zero modulo p. First use 

E Vk = ^rr^(2(5p - 1>' + (p - 2)( P 2 - i)Y), 



fc=i 



24 



then note that for primes p > 3, we have 24\p 2 — 1 and 

12|(p — l)(5p — 1) so that X)fc=o w fc = mod p. Conse- 
quently det(U v ) = 1. 



For the p = 3 case (because of Eq. ( 15 )) we must do a 
little more work to solve a matrix equation analogous to 
Eq. (19). First, we introduce a global phase factor e 1 ^ so 

that 



det (e^£C^ fc2 |fc)(fc| 



1 i.e. 



4-777 
~9~- 



Denote 



a primitive ninth root of unity as ( = e " so that 

det ( c27 ' c (UM]|[,'.])) = 1 

We must permit our qutrit version of [L / 8 to take a more 
general form than that given in Eq. (16), i.e., 

2 

u v = u(v , «!,...)= X! c fc e z 9 ) 

fc=0 

A similar calculation as before leads to the general solu- 
tion (compare with Eq. ( |24[ )) 

(27) 



u = (0, 6z' + 2i + 3e', 6z' + 7' + 6e') mod 9 
For example, letting z' = 1, 7' = 2 and e' = 
w = (wo,fi,W2) = (0, 1,8) 

10 o\ / 1 

=>U V (0,1,8)= I ( = eT 
C 8 / V e- 2 -? 



(28) 



(29) 



One can easily check that all 27 solutions for z', 7', e' € Z p 
obey X)fc=o u fc = mod 3. However, 

det(l^) =££2=0"* = w ( z '+t') 

so that C~ (z ' +y) t/ t , € SU(3). 

We finish this section by noting that knowledge of 
U v D(i\o)U v , U v D( \i-)U v and Eq. ( |Tl"| ) is sufficient to see 
that the effect (modulo an overall phase) of conjugating 
an arbitrary Pauli operator by U v is 



U v Dr x \ z )Ul OC C/r 1 oiir a iv 

Ur/ lJ|Lo:( Z '+2- 1 7'(^-l))+2JJ 



(30) 



C. Group Structure 

For p > 3 the set {U v } C SU(p) with matrix multipli- 
cation satisfies all the prerequisites to be a group: 



Group 


No. 


elements of order 


Min. no. of 


name 


1 


P 


P 2 


P 3 


generators 


p = 2 Z 8 


1 


1 


2 


4 


1 


p = 3 Z 9 x Z 3 


1 


8 


18 





2 


p > 3 


lp 3 




1 





3 



TABLE I. Group structure of the set of diagonal unitaries 
{U v } under matrix multiplication. 



1. Closure: 

U v (zi , 71 , €1 ) U v (z 2 , 72 , £2 ) 
= +z 2 ,7i +72-ei +(-2) 

2. Associativity: Obvious 

3. Identity element: I = 1^(0,0,0) 

4. Inverse element: U^ 1 = U- v 

The fundamental theorem of finite Abelian groups states 
that a finite Abelian group is isomorphic to a direct prod- 
uct of cyclic groups of prime-power order. Furthermore, 
two finite Abelian groups G, G' are isomorphic, G = G' , 
if and only if they have identical structure order i.e., if G 
and G' have the same number of elements of each order. 
We can use this fact to classify the groups that are gen- 
erated by all diagonal gates {U v \i z',-y',e' € Z p } under 
matrix multiplication (see Table [I] for a summary) . For 
p = 2 the gate U^/g is sufficient to generate the entire 
group isomorphic to Zg. For p>3we have 



({Uv}, 



(2t+) 



(31) 



which tells us, amongst other things, that the minimal 
number of generators required to generate this group is 
3. 

For p = 3 one can check explicitly that the gates 
U(z' , 7', e') form a group, and that 

U v (z 1 ,j 1 ,e 1 )U v (z 2 ,j 2 ^2) = (32) 

\U v (zi +Z2,7i +72,^1 +£2) if7i+7 2 <3 
\U v (z 1 + z 2 ,7i + 72, ei +e 2 - 1) if 71 + 72 > 3 

The orders of the individual group elements are 

ord(J7„(0,0,0)) = 1 (33) 
ord(U v (z', 0, e')) = 3 (excluding case @) (34) 
ord(U v (z' , 7', e')) = 9 (excluding cases (33) and |34| ) 

In fact, for p = 3 we have 

({U v }, ■) = Z 9 X Z 3 (35) 

i.e., the group is a direct product of cyclic groups of order 
9 and 3. This group requires only 2 generators. 

The result Eq. ( 15 ) has led to an unusual group struc- 



ture for p = 3 compared to other primes. Combined with 
observations from other authors [r9Tf2"T| on SIC-POVMs 
in p = 3, perhaps it could be argued that, in the context 
of quantum information, 3 is the second oddest prime of 
all 1221. 
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III. GEOMETRICAL FEATURES 



In this section we outline various geometrical proper- 
ties of the gates U v and associated s tates \ ipu„) (as de- 
fined below in Eq. (37l). Subsection III A provides the 
various definitions required for this section. In MB we 



show that the states \ipu v ) are eigenvectors of Clifford 
gates. This may be of independent interest. In |III C| we 
relate both U v and \ipu v ) to the convex hulls of Clifford 
gates and stabilizer states respectively. The reason we do 
so is twofold: (i) We will see how U v and \ipu v ) are singled 
out as being maximally non-Clifford or non-stabilizer in 
some sense, which is interesting in and of itself (ii) The 
preceding geometrical property implies that these gates 
are optimal in the sense of robustness to depolarizing and 
phase damping noise. This latter property should prove 
useful in the context of fault-tolerant quantum computa- 
tion as we discuss in Sec. IIVI 



Useful definitions 



As discussed in Sec. II C the set of gates {U v } forms a 
finite group. In fact, they form a finite subgroup, {U v } C 
{Ug}, of the group of diagonal gates {Ug} defined as 



(0 k e M) 



(36) 



fc=0 



We will often have reason to refer to a state \ipu e ) € C 
that is very naturally associated with the gate Ug via 



Vp 



diag(£/ e ) = Ug 



(37) 



where |+) = (1, 1, . . . , l)J-Jp. A Jamiolkowski state, 

\Ju) € C p2 , corresponding to a unitary operation U G 
U(p) is denoted 



p-i 



\J U ) = (I®U)J2 



m 

o VP 



(38) 




(a) 




J C (Fk\Xk)> 



FIG. 1. Geometry of States and Operations: Filled circles 
(vertices) correspond to stabilizer states or Clifford gates, re- 
spectively, (a) Polytope of Stabilizer States: A single qudit 
has p(p + 1) stabilizer states, the convex hull of which com- 
prises a polytope with p p+1 faces. The states \4>u v } are the 
farthest outside the polytope of all states \ipug) (defined in 
Eq. 371. For p — 2 the stabilizer polytope corresponds to 
an octahedron (as depicted) with 6 vertices and 8 faces. In 
this dimension \ipu v ) = \H) - the magic state introduced 
bv |14l fT5j . Recent results by Campbell et al. [17] show that 
\ipu v ) is a magic state for all p. (b) Polytope of Clifford 
Gates: A schematic picture (the object is actually (p 2 — l) 2 - 
dimensional) of the polytope whose vertices correspond to 
Clifford gates. It is known with certainty in p — 2, and with a 
high degree of plausibility in p £ {3, 5, 7} |16] . that the gates 
U v are farthest outside the Clifford polytope of all U £ U(p). 
Here, each unitary gate, U v or C(f|x) £ C (refer to Eq. (18b) is 
represented by its Jamiolkowski state (as defined in Eq. (|38[)). 



A quantum operation, £ , is a superoperator acting upon 
density operators (i.e., generic quantum states p S T-L p ) 
via 



£ : p in p out i.e., p a 



£ (Pv 



(39) 



The well-known Jamiolkowski isomorphism tells us that 
a complete description of £ is encapsulated in a higher- 
dimensional state Qg defined as 



B. Eigenvectors of Clifford Gates 



As discussed in Sec. |IIB we have an explicit form for 



p" diagonal gates U v (z', 7', e') ( with z',7',e' g Z p ), of 
which p 2 (p—l), corresponding to 7' ^ 0, are non-Clifford. 
Here we prove that each associated state \i>u(z',-y',e')) as 
defined in Eq. ( 37 ) is an eigenvector of the Clifford oper- 
ator 



p-i 



ee = [X®£] [ Y, 



\3,j)(k,k\ 



p 



(40) 



which is the most general form of the operation-state 
duality given in Eq. ( 38 1 . 



C, 



with eigenvalue a/ . Our intuition was that knowledge of 
such eigenstates should prove useful since both |T) and 
I if) (qubit magic states as defined in jTJj) are Clifford 
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eigenstates. The recent result by Campbell et al. [T7] 
whereby states that we call \ipu v ) are shown to be qudit 
magic states, confirms the correctness of this intuition. In 
addition, Zauner's conjecture [33] states that fiducial vec- 
tors of a Weyl-Heisenberg-covariant SIC-POVM lie in the 
eigenspace of a particular class of Clifford gates. While 
the results obtained here are not directly applicable to 
the resolution of the SIC-POVM problem, they may still 
prove useful in this context. 



To prove the claim, first use Eq. ( 18 ) to perform the 
following substitution 

C ( [l 0]|[ i^\ipu(z>,y,e')) = (41) 

^- e U v D {m U^ u(z , tY , e0 ) (42) 

and subsequently use 

Ul\il> u J) = Ul{U v \+)) = \+) 

where |+) is the +1 eigenstate of -D(i|o); by definition. 
Clearly, it follows that 



C 



^ i Oj |[" i Y^uiz'n'^')) = u e \ipu{z',Y,e'))- (43) 





e* D (U v ) 


£*pd(U v ) 


N(\4>u v )) 




p = 


2 45.32% 


14.65% 


0.1036 


2(0.1036)=0.2071 


p = 


3 78.63% 


36.73% 


0.1363 


3(0.1363)=0.4089 


p = 


5 95.24% 


64.00% 


0.1600 


5(0.1600)=0.8000 


p = 


7 97.63% 


73.27% 


0.1202 


7(0. 1202)=0. 8411 



TABLE II. Robustness and negativity: Robustness to noise 
(ej) for depolarizing, e*p D for phase damping) of a gate U 
is the noise rate at which a noisy implementation of U en- 
ters the Clifford polytope. Negativity can be used as a proxy 
for distance outside the relevant (st abi lizer or Clifford) poly- 
tope and is formally defined in Eq. ( 46 1 (states) and Eq. ( 52 1 
(gates). A priori, there is no obvious reason why N(\i/ju v }) 
and N(\Ju v )) should obey such a simple relationship with one 
another. In 116] it was shown that U v were the most robust to 
depolarizing noise of all U £ U(p) (in dimensions 2 to 7 and 
with some caveats regarding an incomplete facet description 
of the Clifford polytope). Here we show (for p £ {2, 3, 5, 7}) 
that U v are also the most robust to phase damping noise of all 
U$. The discussion in Sec. |III C 3| shows that this must also 
imply that \ipUv) are the most robust states (to depolarizing) 
of all states \ipu e )- 



set of all probabilistic mixtures of stabilizer states is an 
interesting geometric object. It is defined as 



C. Noise thresholds for quantum computation 
using U v gates and stabilizer operations 

In this section we show how the gates U v and states 
\ipuS) are exceptional with respect to their convex- 
geometrical relationship to Clifford gates and stabilizer 
states respectively. We will need to define the stabilizer 
polytope (the convex hull of stabilizer states), the Clif- 
ford polytope (the convex hull of Clifford gates) and a 
quantity we call negativity (which can be interpreted as 
a measure of distance outside one of these polytopes). In 
the present context, a state or gate that is farther out- 
side a polytope generally requires more noise (a higher 
degree of impurity) to enter said polytope. To this end 
we introduce a quantity called robustness, which mea- 
sures the amount of noise that can be tolerated before a 
gate (state) becomes expressible as a mixture of Clifford 
gates (stabilizer states). Table [TT| summarizes most of the 
results in this section. How these results were obtained 
is explored in the remainder of this subsection. 



1. Stabilizer Polytope and Clifford Polytope 

For a single-qudit system there are exactly p(p + 1) 
distinct eigenstates of Pauli operators. For p = 2 these 
eigenstates correspond to the vertices of the octahedron 
depicted in Fig.[lja). The Gottesman-Knill theorem tells 
us that a supply of stabilizer states, or mixtures of sta- 
bilizer states, are useless for the task of promoting a sta- 
bilizer circuit to a circuit capable of UQC. As such, the 



p(p+i) 



with < 9i < 1, 9i = 1 



(44) 



where qi can be understood as probabilities and \ip^1 AB ) 
are the aforementioned Pauli eigenstates (stabilizer 
states). 

The convex hull of a finite set of points forms what 
is generally known as a polytope (a higher-dimensional 
generalization of a polyhedron). Using the Minkowski- 
Weyl theorem, we know that every such polytope has a 
description in terms of a finite number of facets (bound- 
ing inequalities, also known as halfspaces). In the present 
context, it is known [21] that exactly p p+1 facets (which 
we denote A) are required to describe STAB. More pre- 
cisely, testing a state p for membership of the polytope 
STAB, leads to the following condition 



p e STAB ^> min Tr [A(u)p] > 0. (45) 



If a state p is outside STAB then we define the negativity 
[27] of the state, N(p) as 



N(p) = \ min Tr [A(u)p}\. 



(46) 



The threshold depolarizing rate, e* D ( 1 4>u e ) ) > 01 a state, 
\ipu e )i is the minimum value of Ed required to make \ipu g ) 
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an element of STAB: 

e* D {\tpv e )) = min£_D (0 < e D < 1) such that (47) 
{l-e D )\i/> Ut )(1> Ua \ + e D ±eSTAB. 

The quantities N(p) and e|j(p) are (inversely) related, as 
discussed in [15] . 

An explicit definition for individual facets A(u) is given 

by 



A(u) = 



1 
P 



n 



(0|l)[«o] 



Xl n (ili-i)K] 



(48) 



where II( |i,)[fc] is the projector onto the ui k eigenspace of 
X a Z b i.e., 



n 



(a\b)[k] 



i (/ + cj- fe x a z 6 + . . . + ^-^(A^zy- 1 ) 

(49) 



Using the Jamiolkowski isomorphisms of Eq. ( 38 ) and 



Eq. (40) we can construct an object (polytope) that is 



analogous to STAB, but where the vertices now corre- 
spond to Clifford gates rather than stabilizer states. As 
before, quantum operations that are expressible as a mix- 
ture of Clifford operations are useless for the task of pro- 
moting a stabilizer circuit to a circuit that is capable 
of UQC. We denote this so-called Clifford polytope as 
CCXTT ps H5] 



CCTTT 



Qs 



j=p( P 2 -i),k=p 2 

Q£= 1i,k\ J c, 
j=i,fc=i 

j =p (p 2 -i).k=p 2 

with < qj <k < 1, ^ Qj,k 
j=i,fc=i 



(^lx fc )^ Jc, (^l«-)' 



(50) 



Testing an arbitrary quantum operation £ for member- 
ship of the Clifford polytope requires construction of 
the associated Jamiolkowski state qs (as described in 
Eq. (40)) and then using 



qs G CCITT 



Tt(Wq £ )>0 (VVKeW), (51) 



where W is a finite set of facets describing CCITF. Anal- 
ogously to Eq. ( 46 ) , we define the negativity of an oper- 
ation 8 as 



N{ ee ) = \ mm Tt[Wq £ ] 



(52) 



The threshold depolarizing rate, e* D {Ue), of a gate, Ug, is 
the minimum value of required to make Ug an element 
of CCITT: 



e* D {Ue) = mmeo (0 < Ed < 1) such that 
{l-e D )\J Ue ){J Ue \+e D ^ G CCXTT. 



(53) 



While W is known to exist, and be finite, the complexity 
of halfspace enumeration is such that we can only claim 
to have derived in [T5] (at least) a subset of W. Never- 
theless, if a given qs (encoding an operation £ ) satisfies 
Tr (W Qs) < for some W G W, then this operation is 
unambiguously outside the Clifford polytope. 



2. Robustness to Depolarizing Noise i.e., 
Maximally Non-Clifford Gates 

In [16] a gate U opt G U(p) was found, for each of 
p G {2, 3, 5, 7}, which required very high amounts of de- 
polarizing noise to become expressible as a mixture of 
Clifford gates. There, it was suggested that the simple 
form of U op t and their high robustness to noise (i.e., high 



e* D in Eq. (58)) made them analogous to the qubit U w /g 
gate in some sense. Here we strengthen the analogy by 
showing that U op t are actually equivalent (i.e., the same 
up to a factor of a Clifford gate) to the gates U v that we 
have derived by enforcing that they should be diagonal 
elements of C3. 

The state \Ju opt ) that is farthest outside CCITT (i.e. 
the convex polytope whose vertices are Clifford gates) is 
that state which achieves 



mm TriWlJuXJul) 
W£W,uesu(p) 



(54) 



where W is the bounding set of facets that describes 
CCXTT. Here we give the explicit relationship between 
highly (and maybe maximally) robust gates given in |16j 
and the generalized versions of C/^/g that we have de- 
scribed in Sec. Ill Bl 



p = 2 
p = 3 



Uopt 

u 



opt = U v C,^-i oj |[o]) ( U v m Eq- VW 



P 
P 



u 



opt 



u v c n _ 



([=i-°i]|[§]) 



(U v in Eq. (26)) 



u 



with U v = U v (z' = 1,7' = 2, e' = 0)] 



The Clifford gates that relate U op t and U v can be ab- 
sorbed into W, creating another W G W with the same 
spectrum, and 



Tr(W'\J Uv )(J Uv \) = Tr (W\J UvC )(Ju v c\) 



(55) 



so that the negativity and robustness results from |16j 
apply here too. 



3. Phase Damping Thresholds via Simplified Jamiolkowski 
Isomorphism 

Phase damping is a physically well-motivated noise 
process (often interpreted as resulting from so-called 



phase kicks), whose overall effect on a state p is to uni- 
formly decrease the amplitude of all off-diagonal elements 
in p (see e.g., 33 ). The implementation of a diagonal 
gate Ug , while suffering phase damping noise (with noise 
rate Epd), results in an overall operation 



£(p) = (1 - e PD )U eP Ul + J2 ( ZkU <>) P i 2 " 



: Ug 



k=l 



= (1 - e PD )U e pUl + ^-j (l - UepUfj (57) 

so that the robustness to phase damping noise, s PD (Ug), 
of a gate, Ug, is 

e PD (Ug) = mm Epd (0 < £ P d < 1) such that (58) 
g £ G CC1FF £ as in Eq. (pi) 



Given a diagonal gate £/# (as defined in Eq. (36l) it is 
trivial to see that the state \ipu g ) ( as defined in Eq. (37)) 
provides a complete description of the gate. A quan- 
tum process consisting of a probabilistic mixture of var- 
ious different diagonal gates Ug is thus representable by 
a quantum state 



0<qi<l,J29i 



(59) 



Effectively, this is a simplified form of the Jamiolkowski 
isomorphism (given in Eq. (40)) that is only possible be- 



cause our allowed operations are highly restricted. With 
that said, the operation in Eq. (56) can equally well be 



represented by a single-qudit state: 



P = (1 - £PD)Hu e )(lpUe\ + 51 I^Z fc t/ e )(V'Z 



k=l 



= (i - e PD MuMu a \ + —Ai- \i>uMu e \) 

p- 1 

To summarize: a noisy (phase-damped) implementation 
of Ug can be identified as a state of the form Eq. (59). 
As discussed in Sec. II B there are exactly p 2 diago- 



nal Clifford gates, corresponding to U v (z' , 0, e'), and each 
such gate corresponds to a state \4>u v )- These p 2 states 
comprise the vertices of a polytope of dimension p(p — 1) 
contained in the space spanned by Eq. (59). For exam- 
ple (and consulting Fig. [lja)), in the p — 2 case one 
sees that the (two-dimensional) x-y plane of the octahe- 
dron contains p 2 — 4 vertices. Checking the negativity 
(distance outside this polytope) of a given state amounts 
to evaluating its expectation value with respect to all 
p p — 4 facets that comprise the polytope boundary (i.e., 
the four edges of the octahedron that are contained in the 
x-y plane). Of all states \tpu e ), the one farthest outside 
this polytope is iV'E/^/g)- The general expression for all 
p p distinct facets of this p(p — l)-dimensional polytope 



(which we call ^4 e d ge ) is 



p-i 



A edge (u 1 ,u 2 , ...,u p ) = - ^2 A(u ,ui,. ..,u p ) (60) 



u =0 



'4- 



> i i ' »,] - - I ( 61 ) 



(56) where A are the facets of STAB as defined in Eq. (48) 



The final element that is required is to realize that an 
operation Ug that is maximally robust to phase damping 
noise is exactly the operation for which | ipjj e ) is most re- 
sistant to depolarizing noise before entering the p(p — 1)- 
dimensional polytope discussed above. In fact, a simple 
calculation shows that 



-PD 



{Ug 



P-I 
P 



where e PD (Ug) is the phase-damping noise rate required 
to make Ug enter the convex hull of diagonal Clifford 
gates, and SptWug}) is the depolarizing noise rate re- 
quired to make \4>u e ) enter the convex hull of stabilizer 
states. 

Any facet ^4 d g e can be decomposed as 



A edge = ^A i |A j XA,-| (Ai<A 2 ...), 



(62) 



and the state |Ai) is the state that, out of all qudit states 
p, maximally violates Tr [A G d gc p] > 0. A simple calcula- 
tion shows that the minimum eigenvalue of ^4 d g c (over 
all possible A c d gc ) is Ai(^4 c d go ) = — (p — l)/p 2 > when p is 
an odd prime. Consequently, if Ai = — (p — l)/p 2 , and 
|Ai) is of the form \ipug), then |Ai) is maximally robust to 
depolarizing noise and the operation Ug that |Ai) repre- 
sents is maximally robust to phase-damping noise. This 
is the case for p € {2, 5} in our current investigation. For 
p G {3, 7}, the states |Ai) that achieve Ai = — (p — l)/p 2 
are not of the form \ipu e ) aud so we had to resort to a 
numerical optimization over all ^4 e d g e aud all states \ipUe)- 
For p = 3 there are two distinct types of edge, as clas- 
sified by spectrum: 

A(A:dgc) = {— g, g, 9} 



or 



A(A cdgc ) = I ~ (3 sin i - V3 cos g 
- + 3 sin - cos ~ 
i(l + 2V3cos^)} 



There are p 2 (p — 1) = 18 of the latter A cdgc , where 
the minimizing eigenvector for each distinct facet cor- 
responds to a distinct non-Clifford U v . In p = 5 there 
are exactly p 2 (p — 1) = 100 edges with spectrum 

A(A cdgc ) = {-0.16,-0.08361,0.04,0.04,0.36361} 
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and these correspond to the 100 non-Clifford U v . Inp = 7 
there are at least 2(7 2 ) = 98 facets A c a ge with minimal 
eigenvalue —0.12016, whose corresponding eigenvector is 
of the form \il>u v )- These were the most robust states 
of all \ipu a ) that we could find by optimization, but it is 
possible that we became trapped in a local minimum. 

As a final comment, we note that qubit form of this ar- 
gument (simplified Jamiolkowski isomorphism etc.) was 
presented by Virmani et al. [29] (see also [30]) where an 
adversarial phase damping model was used to obtain up- 
per bounds on the quantum fault-tolerance threshold. In 
that case it was found that 



the form 



[14.7% « £* PD (U w/8 )] 



(\ifo w , a )) * 29.3%] 



(63) 



which implies that the U w /g gate, whilst maximally ro- 
bust amongst all diagonal gates, requires about 15% 
phase damping noise before it becomes expressible as a 
mixture of (diagonal) Clifford gates. 



p-i 

3=0 



with A, = p 



p + 1 

4 



Note that 



M (3) = U v (z' 
M^ = U v (z' 
MW = u v (z' 



l, 7 ' = l,e'=0) 

2, y = l,e' = 2) 

3, y = l,e / =4) 

rfr)\ 



(67) 
(68) 



(69) 
(70) 
(71) 



in our notation. They showed that \Mq ), which corre- 
spond to \i[>Uv) defined in Eq. (37), are distillable for all 
prime dimensions and hence are magic states. For p = 3 
a distillation routine with remarkably good performance 
was found and we use this result in the discussion below. 
If we have access to a superoperator £m u \ e ) defined 

as 



IV. APPLICATIONS IN FAULT-TOLERANT 
QUANTUM COMPUTING 

In [T7] Campbell et al. quote results by Nebe, Rains 
and Sloane j34J [35] which show that the gate set 



(C, CSUM, [/) 



(64) 



is dense in SU(p"), where C is the set of single-qudit 
Cliffords, CSUM is the generalized version of the CNOT 
gate and U is a non-Clifford single-qudit gate, 



CSUM : |o)|6) h> |o)|a- 
U e SU(p)\C. 



b mod p) 



(65) 
(66) 



In particular, U v is sufficient to promote multi-qudit Clif- 
ford gates to a universal gate set. We expect that U v 
possesses all the additional qualities that makes U w / S the 
preferred non-Clifford gate in qubit-based universal gate 
sets. 

In [5] [TS] it is argued, for the qubit case, that cre- 
ation of a Clifford eigenstate should be easier to do in 
a fault-tolerant manner than fault-tolerant implementa- 
tion of U v directly. It is easy to see that, for an arbitrary 
qudit state \ip ai b), 

n (0 ,0|l,d-l)[0](l^(7j ® IV'arb)) = (^l^arb)) ® |0) 

where n( 0i o|i,d-i)[o] denotes a rank-p projector onto the 
lo° eigenspace of the operator Z <£> Z~ x (i.e., a stabilizer 
measurement). Clearly, creation of Clifford eigenstates 
\ipuj) (see Sec. Ill B I is sufficient to promote a stabilizer 



circuit to UQC in the qudit case too. 

In [T7] Campbell et al. introduce a qudit gate M of 



£ Wv v ),e)(j>) = (l-e)U v pUl+i 



(72) 



then it can be used to create noisy versions of \ipu v ) m 
a straightforward way, as depicted in Fig. 2^a). How- 
ever, a simple circuit given in Fig. [2jb) shows how a 
less noisy version of \ipu v ) can be created by using the 
same operation }, e ) as well as some additional sta- 

bilizer operations. The straightforward method produces 
a state \if>u v ) with effective depolarization rate e whereas 
the postselected version produces \i^uj) with effective de- 
polarization rate e'{< e). The relationship between the 
effective noise rates is 



pe 



and, inverted, 



1 + (p - l)e> 



p- (p- l)e' 



(73) 



(74) 



In |17| it was shown that, for p = 3, depolarized ver- 
sions of \ipu v ) could be distilled for noise rates up to about 
32%. The implication is that a superoperator ), e ) 
enables universal quantum computation (via MSD) up to 
noise rates of around 58%. This is found by solving 



Our results in Sec. 
can never enable 



IIIC2 



0.3165 
0.5815 



and 1161 indicate that £ , 



(75) 
(76) 

<Wtr„>,<0 



universal quantum computation (by 
supplementing Clifford gates) for noise rates e > 78.6%. 
In the qubit case it was shown, using similar argu- 
ments, by Reichardt[2I] and Buhrman et al. [3S] (and 
more generally in 32J), that the noise rates for which 
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1+) 



- 0--e)WuMu v \+e; 

(a) 



1+) 



l°>-© 



n 



(b) 



10) 



(l-e')\^ Uv ){^ Uv \+e'- 



FIG. 2. Dilution of Noise in Magic State Preparation: (a)Straightforward Magic State Preparation: Using the superoperator 
£(\i>u ),e) to prepare an imperfect (depolarized) version of a magic state \ipu v )- (b)Postselected Magic State Preparation: Using 
the same superoperator £n^, v ), e ) to create a less imperfect (e' < e) version of the magic state \ipu v }- The Pauli measurement 
operator II stands for postselection on receiving the outcome +1 using the measurement II(c lom, p -i)[o] (see Eq. (67l). The 
circuit elements to the left of II implement the creation of the Jamiolkowski state ge (see Eq. (40 1) describing £{\^ u >, e ). 



Lower Bound Upper Bound 



in Tabic. |III| can potentially be lowered. 



p = 


2 


45.32% 


45.32% 


p = 


3 


58.15% 


78.63% 


p = 


5 


80.61% 


95.20% 


p = 


7 


72.24% 


97.63% 



TABLE III. Noise thresholds for universal quantum compu- 
tation: When supplementing the set of Clifford gates, which 
are presu med to be perfect, with a depolarized version of U v 
(as in Eq. |72|), then there exist regimes of the noise param- 



eter, s, for which UQC is either provably possible or provable 
impossible. Noise rates between these two bounds are those 
for which we currently have no proof regarding the utility 
of depolarized U v when supplementing stabilizer operations. 
Lower bounds are evaluated by using the results of Camp- 
bell et al. 1 17j . along with the noise dilution protocol of Fig. [5] 
and Eq. (74 1. Upper bounds are established via geometrical 
arguments in Sec. Ill C 2 and 16 . 



,>.<0 



was outside the Clifford polytope were ex- 



actly those noise rates for which £(\^ v g > >£ ) could sup- 
plement Clifford gates to enable universality. It is an 
interesting open question for qutrit systems whether 
the operation £(U V ), s ) with noise rates in the range 
58.15% < e < 78.6% can enable universal quantum com- 
putation. All of the above techniques can, in principle, be 
applied to any prime dimensional \ipu v ) and we present 
a summary in Table. Ill for dimensions p € {2,3,5,7}. 
Indeed, the performance of the p = 5 MSD routine is 
such that the gap between upper and lower bounds in 
Table. |III| is even smaller in this case. If additional qu- 
dit MSD routines are developed, then the lower bounds 
of Table. |III| can potentially be raised. Similarly, if an 
analagous scenario to that of bound non-stabilizer states 
[2"6"ll3"rj] also holds for operations, then the upper bounds 



V. SUMMARY & OPEN QUESTIONS 



Motivated by the utility and geometric prominence of 
the qubit U^/s gate, we provided an explicit solution 
for all diagonal qudit gates that displayed the same re- 
lationship with the Clifford group (i.e., we constructed 
diagonal gates from the third level of the Clifford hi- 
erarchy). We saw that these diagonal gates generated 
a finite group whose structure depended upon whether 
p = 2, p = 3 or p > 3. It might be interesting to fully 
enumerate all the single-qudit elements of C3, or analyze 
the structure of the diagonal subset of Ck (particularly 
for p — 3). Geometrically, these generalized U n /g gates, 
which we have called U v , appear to display the same re- 
lationship with the set of Clifford gates - a relationship 
which makes them maximally non-Clifford in some sense. 
The state \tpu v ) S C p , defined as = U v \+), was al- 

ready known to be useful and geometrically significant in 
the p = 2 case (where it is widely known as the | if) -type 
magic state), and we discussed some properties of the 
general qudit case which led us to believe they could also 
be useful. As we completed this work we became aware 
of results by Campbell et al. [T7] which show that states 
\tpu v ) a re indeed magic states for all prime dimensions. 
In the final section we use the results of Campbell et al. 
to show noise rates for which noisy versions of U v can 
and cannot provide UQC (when supplementing the full 
set of Clifford gates). A very interesting problem is to 
further close the gap between these noise regimes, a gap 
that is non-existent in the qubit case. 
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